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We study the properties of error correcting codes for noise models in the presence of asymmetries 
and/or correlations by means of the entanglement fidelity and the code entropy. First, we consider 
a dephasing Markovian memory channel and characterize the performance of both a repetition code 
and an error avoiding code {Crc and Cdfs, respectively) in terms of the entanglement fidelity. 
We also consider the concatenation of such codes {Cdfs ° Crc) and show that it is especially 
advantageous in the regime of partial correlations. Finally, we characterize the effectiveness of the 
' codes Cdfs, Crc and Cdfs ° Crc by means of the code entropy and find, in particular, that the 

(N : effort required for recovering such codes decreases when the error probability decreases and the 

pH ' memory parameter increases. Second, we consider both symmetric and asymmetric depolarizing 

I noisy quantum memory channels and perform quantum error correction via the five qubit stabilizer 

code C[[5_i_3]]. We characterize this code by means of the entanglement fidelity and the code entropy 
' as function of the asymmetric error probabilities and the degree of memory. Specifically, we uncover 

that while the asymmetry in the depolarizing errors does not affect the entanglement fidelity of 
the five qubit code, it becomes a relevant feature when the code entropy is used as a performance 
quantifier. 

PACS numbers: Quantum Error Correction (03.67.Pp); Decoherence (03.65. Yz). 
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I. INTRODUCTION 



Quantum error correction (QEC) is a theoretical scheme developed in quantum computing to defend quantum 
coherence against environmental noise. There are different methods for preserving quantum coherence. One possible 
technique exploits the redundancy in encoding information. This scheme is known as "quantum error correcting codes" 
(QECCs). For a comprehensive introduction to QECCs, we refer to Within such scheme, information is encoded 
] in linear subspaces (codes) of the total Hilbert space in such a way that errors induced by the interaction with the 
' environment can be detected and corrected. The QECC approach may be interpreted as an active stabilization of a 
quantum state in which, by monitoring the system and conditionally carrying on suitable operations, one prevents the 
loss of information. In detail, the errors occur on a qubit when its evolution differs from the ideal one. This happens 
by interaction of the qubit with an environment. Another possible approach is the so-called noiseless codes (also 
known as decoherence free subspaces (DFSs) or error avoiding codes). For a comprehensive introduction to DFSs, we 
refer to 0]. It turns out that for specific open quantum systems (noise models in which all qubits can be considered 
symmetrically coupled with the same environment), it is possible to design states that are hardly corrupted rather 
than states that can be easily corrected (in this sense, DFSs are complementary to QECCs). In other words, it is 
. , possible to encode information in linear subspaces such that the dynamics restricted to such subspaces is unitary. This 
. implies that no information is lost and quantum coherence is maintained. DFS is an example of passive stabilization 
■ ■ ' of quantum information. 

The formal mathematical description of the qubit-environment interaction is usually given in terms of quantum 
channels. Quantum error correction is usually developed under the assumption of i.i.d. (identically and independently 
distributed) errors. These error models are characterized by memoryless channels A such that n-channel uses is given 
by A(") = A®". In such cases of complete independent decoherence, qubits interact with their own environments 
which do not interact with each other. However, in many physical situations, qubits do interact with a common 
environment which unavoidably introduces correlations in the noise. For instance, there are situations where qubits 
in a ion trap set-up are collectively coupled to their vibrational modes In other situations, different qubits in a 
quantum dot design are coupled to the same lattice, thus interacting with a common thermal bath of phonons 
The exchange of bosons between qubits causes spatial and temporal correlations that violate the condition of error 
independence Q. Memory effects then arise among channel uses with the consequence that A(") ^ A«5". Recent 
studies try to characterize the effect of correlations on the performance of QECCs |6l-[ll|. It appears that correlations 
may have negative or positive Q impact on QECCs depending on the features of the error model being considered. 

A part from being correlated, noise errors may also be asymmetric. Most of the quantum computing devices [T^ 
are characterized by relaxation times (rreiaxation) that are one- two orders of magnitude larger than the corresponding 
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dephasing times (rdcphasing)- Relaxation leads to both bit- flip and phase- flip errors, whereas dephasing (loss of 
phase coherence) only leads to phase-flip errors. Such asymmetry between Ti-daxation and Tdophasing translates to an 
asymmetry in the occurrence probability of bit- flip (px) and phase-flip errors (pz)- The ratio ^ is known as the 
channel asymmetry. Quantum error correction schemes should be designed in such a way that no resources (time and 
qubits) are wasted in attempting to detect and correct errors that may be relatively unlikely to occur. Quantum codes 
should be designed in order to exploit this asymmetry and provide better performance by neglecting the correction 
of less probable errors [l3l - [l5| . Indeed, examples of efficient quantum error-correcting codes taking advantage of 
this asymmetry arc given by families of codes of the Calderbank-Shor-Steane (CSS) type (asymmetric stabilizer CSS 

codes) 

Devising new good quantum codes for independent, correlated and asymmetric noise models is a highly non trivial 
problem. However, it is usually possible to manipulate in a smart way existing codes to construct new ones suitable 
for more general error models and with higher performances [isl IToj . Concatenation is perhaps one of the most 
successful quantum coding tricks employed to produce new codes from old ones. The concept of concatenated codes 
was first introduced in classical error correcting schemes by Forney [20j . Roughly speaking, concatenation is a method 
of combining two codes (an inner and an outer code) to form a larger code. In classical error correction, Forney has 
carried on extensive studies of concatenated codes, how to choose the inner and outer codes, and what error probability 
threshold values can be achieved. The first applications of concatenated codes in quantum error correction appear in 
pll [22} . In the quantum setting, concatenated codes play a key role in fault tolerant quantum computation and in 
constructing good degenerate quantum error correcting codes. 

Obviously enough, the choice of a code depends on the error model. However, in selecting a code rather than 
another, one should take into account their performance on the error model. Since there is no single quantity holding 
all information on a code, different quantum code performance quantifiers have appeared into the QEC literature 
(23j . For instance, in [23| it is suggested that the entanglement fidelity [2^ is the important quantity to maximize 
in schemes for quantum error correction. In j26j . it is proposed that the code entropy is an important quantitative 
measure of the quantum correction schemes by determining a hierarchical structure in the set of protected spaces [l^j . 

Motivated by all the above mentioned considerations, in this work we study the properties of error correcting codes 
for noise models in the presence of asymmetries and/or correlations by means of the entanglement fidelity and the 
code entropy. First, we consider a dephasing Markovian memory channel and characterize the performance of both a 
repetition code and an error avoiding code {Crc and Cdfs^ respectively) in terms of the entanglement fidelity. We 
also consider the concatenation of such codes {Cdfs ° Crc) and show that it is especially advantageous in the regime 
of partial correlations. Finally, we characterize the effectiveness of the codes C^ps^ Cue and Cjjps o Crc by means 
of the code entropy. We show that the effort required for recovering such codes decreases when the error probability 
decreases and/or the memory parameter increases. Second, we consider an asymmetric depolarizing noisy quantum 
memory channel and perform quantum error correction via the five qubit stabilizer code C[[5 1 3]] . We characterize 
this code by means of the entanglement fidelity and the code entropy as function of the asymmetric error probabilities 
and the degree of memory. Specifically, we uncover that while the asymmetry in the depolarizing errors does not affect 
the entanglement fidelity of the five qubit code, it becomes a relevant feature when quantifying the code performance 
by means of the code entropy. 

The layout of this article is as follows. In Section II, we present few theoretical tools employed in the rest of 
the article. In particular, we emphasize the conceptual and operational meanings of the entanglement fidelity and 
code entropy, the two code performance quantifiers employed here. In Section HI, we study a dephasing Markovian 
memory channel and characterize the performance of both a repetition code and an error avoiding code {Crc and C^ifs, 
respectively) in terms of the entanglement fidelity. We also consider the concatenation of such codes {Cdfs°Crc) and 
show that it is especially advantageous in the regime of partial correlations. Finally, we characterize the effectiveness 
of the codes Cdfs, Crc and Cdfs°Crc by means of the code entropy. In Section IV, we analyze both symmetric and 
asymmetric depolarizing noisy quantum memory channels and perform quantum error correction via the five qubit 
stabilizer code C[[5 1 3]]. We characterize this code by means of the entanglement fidelity and the code entropy as 
function of the asymmetric error probabilities and the degree of memory. We find, among other things, that while 
asymmetric depolarizing errors do not affect the entanglement fidelity of the five qubit code, they do affect its code 
entropy. A summary of the uncovered results and our comments on the usefulness of both the entanglement fidelity 
and code entropy as suitable code performance quantifiers for the error models considered appear in Section V. 

II. PRELIMINARIES 

In this Section, we briefly describe some basic features of quantum error correcting codes and decoherence free 
subspaces. Finally, we point out the conceptual and operational meanings of the entanglement fidelity and code 
entropy, the two code performance quantifiers employed in this work. 
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A. On quantum error correction schemes 

1. Quantum error correcting codes 
The formal mathematical description of a system-environment interaction is usually given in terms of quantum 



channels [23j. When a quantum system in a state p (belonging to the set of density operators a{'H) defined over the 
Hilbert space T-L of dimension d associated to the system) is exposed to the interaction with the environment, the 
result is the noisy state describable by the action of a completely positive and trace preserving map A : (j{H) — ?> cr{H) 
with, 



Aip) = J2AkpAl. (1) 



The operators : % ^ % (Kraus operators) define the errors affecting the system's state p. It is well known [2j] 
that a standard quantum error correction scheme is characterized by a subspace C C 'H such that for some Hermitian 
and positive matrix of complex scalars F = (7;^.) the corresponding projection operator Pq fulfills the error correction 
condition, 

PcAlA„,Pc^^hnPc, (2) 

for TO = 1,..., (P. The matrix F is also known as the error correction matrix. The action of the correctable error 
operators {Ajconoctabic ^ {^k} on vectors spanning C (codewords) allows to define the (completely positive and 
trace preserving) recovery map TZ ^ {Ri}- The composition of the recovery map TZ with the error map A gives the 
recovered channel, 



Arcc. (p) (7^ ° A) (p) = ^ ^ {RiAk) p (RiAk)^ . (3) 



I k 

For a comprehensive introduction to QECCs, we refer to 



2. Decoherence Free Subspaces 

Following 0, we mention few relevant properties of DFSs. Consider the dynamics of a closed system composed of 
a quantum system Q coupled to a bath B. The unitary evolution of the closed system is described by the combined 
system-bath Hamiltonian i/tot, ^^tot = Hq ^Ib + Hb ®Iq+ -ffint, H-mt = Yla Ea^Ba- The operator Hq {Hb) is the 
system (bath) Hamiltonian, Iq [Ib) is the identity operator of the system (bath), Ea are the error generators acting 
solely on Q while act on the bath. The last term in TJtot is the interaction Hamiltonian. 

A subspace T^dfs of the total system Hilbert space is a decoherence free subspace if and only if: i) Ea\i}) = Cq 
with Cq, € C, for all states \'ip) spanning HDrs, and for every error operator in -ffint- In other words, all basis states 
spanning T-Ldfs a-rc degenerate eigenstates of all the error generators Ea] ii) Q and B are initially decoupled; iii) Hq |^/;) 
has no overlap with states in the subspace orthogonal to 'Hdfs- To establish a direct link between QECCs and DFSs, 
it is more convenient to present an alternative formulation of DFSs in terms of the Kraus operator sum representation. 
Within such description, the evolution of the system Q density matrix is written as, pg (t) = Trg \JJ {pQ (E) pb) = 

J2k^i^PQ (0) where U — e~*^'°'* is the unitary evolution operator for the system-bath closed system and the 
initial bath density matrix pB equals ^„ A„ \n) {n\. The Kraus operators Ak (satisfying the normalization condition) 
are given by Ak = yfKi {m\U\n) with ^l^fc = Iq with k = (n, m) and where \m) and \n) are bath states. It turns 
out that a A^DFS-dimensional subspace T^dps of H is a DFS if and only if all Kraus operators have an identical unitary 
representation (in the basis where the first A^dfs states span HoFs) upon restriction to it, up to a multiplicative 
constant, 

where gk — ^Jn {m\Uc\n) and Uc = g^*^"^* with He = Hb +i?int- Furthermore, Ak is an arbitrary matrix that acts on 

^DFS (with H — Hbfs ffi^iiFs) ^^'^ ^^^y cause decoherence there; Uq'^^^ is Uq restricted to 'Hdfs- Now recall that 
in ordinary QECCs, it is possible to correct the errors induced by a given set of Kraus operators {Ak} if and only if, 

Kkic 

•fe 



RrAk = ( ) , (5) 
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Vr and k, or equivalently, 



where {Rr} are the recovery operators. The first block in the RHS of ([S]) acts on the code space C while the matrices 
Brk act on C""- where H = C © C""-. From (j4|) and ([5]), it follows that DFS can be viewed as a special class of QECCs, 
where upon restriction to the code space C, all recovery operators Rr are proportional to the inverse of the system Q 

evolution operator, Rr oc (^C^g^^^''^ ■ Assuming that Rr is proportional to the dagger of Uq'^^^ , from ^ and © it 

also turns out that A^. cx J7g^^^^ upon restriction to C. Furthermore, from Q and ([5]), it follows that 7^.^./ = g^gk'- 
However, while in the QECCs case jkk' is in general a full-rank matrix (non-degenerate code), in the DFSs case this 
matrix has rank 1. In conclusion, a DFS can be viewed as a special type of QECC, namely a completely degenerate 
quantum error correcting code where upon restriction to the code subspace all recovery operators are proportional to 
the inverse of the system evolution operator. As a side remark, in view of this last observation we point out that it is 
reasonable to quantify the performance of both active and passive QEC schemes by means of the same performance 
measure. In what follows, we will briefly describe the conceptual and operational meanings of the entanglement 
fidelity and code entropy, respectively. 



B. On code performance quantifiers: entanglement fidelity and code entropy 

1. Entanglement Fidelity 

Entanglement fidelity is a useful measure of the efficiency of QECCs. It is a quantity that keeps track of how 
well the state and entanglement of a subsystem of a larger system are stored, without requiring the knowledge of 
the complete state or dynamics of the larger system. More precisely, the entanglement fidelity is defined for a mixed 
state p = "^^iPiPi ^^^Hr IV') (V'l in terms of a purification \ip) S H (g) Hr to a reference system Hr. The purification 
\ip) encodes all of the information in p. Entanglement fidelity is a measure of how well the channel A preserves the 
entanglement of the state H with its reference system Hr. The entanglement fidelity is formally defined as follows 

^ ip, A) "^'{Mi^® Inn) mmii'), (7) 

where \ip) is any purification of p, I-Ur is the identity map on a (1-Lr) and A® is the evolution operator extended 
to the space T-L^T-Lr, space on which p has been purified. If the quantum operation A is written in terms of its Kraus 
operator elements {Ak\ as, A [p] = AkpAl., then it can be shown that the operational expression of (O becomes 



-F (p. A) = ^ tr {A,p) tr (4p) = ^ |tr ipA,)f . (8) 



k k 

This expression for the entanglement fidelity is very useful for explicit calculations. Finally, assuming that 

A:(j{H)3p^A{p)^Yl AkpAl ea{n), diuicH =N (9) 

and choosing a purification described by a maximally entangled unit vector 1-0) S H <^ H for the mixed state p = 
dhkn^^ , we obtain 



J-( A) = 



fe 



^EM'^I'- (10) 



The expression in (jlOp represents the entanglement fidelity when no error correction is performed on the noisy channel 
A in 
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2. Code entropy 

A convenient quantifier of the action of the map A in ([l} on a initial quantum state p is represented by the entropy 
change 5 (A (p)) — S (p). It can be shown that [28j . 

< |5 (A (p)) -S{a)\<S{p)<S{A (p)) + S {a) , (11) 

where a — a (A, p) is the so-called Lindblad matrix, an auxiliary quantum state with matrix elements aim defined as, 

='Tr(p4A„) , (12) 

with I, m — 1, 2,.., (P. The bounds provided in (fTTj) are the so-called Lindblad bounds [28!]. The quantity S {a) is 
known as the entropy exchange of the operation A. From (jlip . it follows that if p is pure, S (p) = and therefore 
Sia) = S{A{p)). 

As pointed out earlier, a standard quantum error correction scheme for a given quantum operation A is characterized 
by a subspace C such that for some Hermitian and positive error correction matrix of complex scalars T = (jim) the 
corresponding projection operator fulfills the relation A quantum error correcting code for A is determined by 
the subspace related to Pc- The rank of F is bounded above by the dynamical Choi matrix 2?a associated with A (p), 

V^^^'Tri^AlAm). (13) 

For orthogonal Kraus operators {Ai}, the Choi matrix equals diSim with < d/ representing the eigenvalues of Da. 
The rank of 2?a equals the minimal number of Kraus operators needed to describe A. Given a code C for A and 
assuming that the initial quantum state p belongs to the code subspace (that is, PcpPc ~ P)i it turns out that 
c/m = lim IH]- In other words, the error correction matrix F equals the Lindblad matrix a (A, p) provided that the 
initial quantum state p belongs to the code subspace. Motivated by these considerations, given a quantum operation 
A with Kraus operators {A{\ and a code C with error correction matrix F, Kribs and coworkers named the von 
Neumann entropy S (A, C), 

5(A,C)''=i'5(F) = -Tr(FlogF), (14) 

the entropy of C relative to A [2^. They show that, 

< 5(A, C) < logZ?, (15) 

where S (A, C) = iff C is a unitarily correctable code for A and S (A, C) = logi? iff C is a non-degenerate code for 
A. From a conceptual point of view, the entropy of a code C relative to a quantum noisy channel A can be viewed 
as a measure quantifying the nearness of the given error correcting code C to a decoherence free subspace Cdfs • The 
closer C is to a Cdfs, the smaller is the code entropy 5 (A, C). The smaller is the code entropy, the smaller is the 
amount of effort required to recover the quantum state corrupted by the noise by means of a suitable error recovery 
operation. The simplest scenario occurs for codes with zero entropy. Such codes are known, as we said, as unitarily 
correctable codes and can be recovered with a single unitary operation. In particular decoherence free subspaces are 
a special class of unitarily correctable codes where the recovery is given by the trivial identity operation. 

We recall that in standard active error correction schemes, the action of the recovery operation pushes all the noise 
into the ancillary qubits, so that errors are eliminated when the ancilla is traced out (29[. That said, we emphasize 
that the numerical value of the code entropy quantifies the number of ancilla qubits needed to perform a recovery 
operation [26| and the rank of F gives the number of Kraus operators necessary to describe the action of A restricted 
to C and therefore the number of Kraus operators necessary for a recovery operation. 

III. MODEL I: DEPHASING MARKOVIAN MEMORY CHANNEL 

In this Section, we consider a dephasing Markovian memory channel and characterize the performance of both a 
repetition code and an error avoiding code {Crc and Cdfs, respectively) in terms of the entanglement fidelity. In 
particular, we consider the concatenation of such codes {Cdfs ° Crc) and show that it is especially advantageous in 
the regime of partial correlations. Finally, we characterize the effectiveness of the codes Cdfs, Crc and Cdfs ° Crc 
by means of the code entropy. We show that the effort required for recovering such codes decreases when the error 
probability decreases and the memory parameter increases. 
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A. Entanglement fidelity-based analysis 

Repetition code for correlated phase flips. The model considered is a depliasing quantum Markovian memory channel 
A''")(p). In explicit terms, we consider n qubits and Markovian correlated errors in a dephasing quantum channel, 



A'"'(p)=^ P^r.\^r.-lP^„-^l\^r.-2■■■P^2\^lP^A^^^®■■■^^^l)pi^^^®■■■^^^l)\ (16) 



dcf 

ii ^7l— 



where Aq =^/, Ai ^= Z are Pauli operators. Furthermore the conditional probabilities Pi^\i^_^ are given by, 

= (1 - /-Op*. ^k^i: P^k=0 = l-p,Pt, = l ^P, (17) 

with, 

1 

X! P^,.K-iP^,^-iK-2---Pt2i.,Pii (18) 

21,..., in—0 

To simplify our notation, we may choose to omit the symbol of tensor product "(8)" in the future. Ai^^...(E}Ai-^ = 
Ai^...Ai-^. Furthermore, we may choose to omit the bar "|" in Pi^\i^ and simply write the conditional probabilities 
as Pi^i . QEC is performed via the three-qubit repetition code. Although the error model considered is not truly 
quantum, we can gain useful insights for extending error correction techniques to fully quantum error models in the 
presence of partial correlations. The performance of quantum error correcting codes is quantified by means of the 
entanglement fidelity as function of the error probability p and degree of memory /x. By considering the case of 
with rt = 3, it follows that the error superoperator A associated to channel is defined in terms of the following error 
operators, 

A ^ {A',.,.., 4} with A(3)(p) "JL' ^ A',pA2 and, ^ A'^A', ^ h^s- (19) 

fe=0 fc=0 

In an explicit way, the error operators {Ag,.., Aj} are given by, 

A', = v^F^' ® ® = vpF^' I^ A', = v^/i ® z2 ® /3, 

A'^ = \f^h^ 0/2® Z^, A'i = ypf^Z^ ®Z^® /3, A'^ = \IW^Z^ 0/2® Z^, 



A'q = \J p'i' r ® Z^ ® Z\ A't = \l p\''' Z^ Z' Z-^ , (20) 

(3) 

where the coefficients for fc = 1,.., 7 read, 

~(3) 2 ~(3) ~(3) ~(3) 

Pa = PooPo, Pi = POOPlOPO, P2 = POlPlOPO, P3 = PooPoiPi, 

~(3) -(3) -(3) ~(3) 2 

P\ = PlOPllPO, Ph = PoiPwPi, Ph = POlPllPl, P) = PiiPl, (21) 

with, 

Po = (1 - P) , Pi = P, Poo = ((1 - m) (1 - P) + m) , 

Poi = (1 - A^) (1 -p) , Pio = (1 - m)p, Pu = ((1 - Ai)P + Ai) ■ (22) 
Then consider a repetition code that encodes 1 logical qubit into 3-physical qubits. The codewords are given by, 

|0)^|0l) = + |1)^|1l>'= I ), (23) 

where |±) =^ ^"^^"'"^ . The set of error operators satisfying the detectability condition [sO], PcA'j^Pc = Xa'^Pc, where 
Pc = |0l) (Oil + |1l) (1l| is the projector operator on the code subspace C = Span{\QL) , |1l)} is given by, 

-Adctcctabic = {A'q, A[, 4, 4, A^, A'e} C A. (24) 



7 



The only non-detectable error is A'j. Furthermore, since all the detectable errors are invcrtible, the set of correctable 
errors is such that ^^Qj-roctabic-^concctabic is detectable. It follows then that, 

-^correctable — {^0' ^1' ^2' — "'^detectable '!= "4. (25) 

The action of the correctable error operators .4correctabic on the codewords |0l) and is given by, 

\0l) a^|Ol) = VpF|+ + +>,^;|Ol> = v^|- + +), A^|OL> = v/^|+-+),^:,|Oi) = VpF|+ + -) 



\H) ^ A', lU) = I ) , A[ - ^i^> 1+ - -) , A', |U) = ^pt> I- + -) , A', = V^r' I- - ^26) 

The two four-dimensional orthogonal subspaccs V"^ and V^^ of H2 generated by the action of .4correctabie on |0l) 
and result, 



Span{\v',-) = \+ + +),\v",-)^\- + +), |«3"") = |+-+), |t;!;-) = |+ + -) }, 



(27) 



and, 



V'- =Span{\vl-) = \ ).\vl^) = \+--), H') = h + -) , H') = h - +)} , (28) 

respectively. Notice that V"'- ® V^^ = 'H2- The recovery superoperator TZ ^ {Ri} with I = 1,..,4 is defined as [3l| . 



(29) 



where the unitary operator Vi is such that V/ = for i £ {0, 1}. Substituting ([27]) and (|28p into (|29p . it 

follows that the four recovery operators {Ri, R2, R3, R4} are given by, 

Ri = \0l) (OlI + lU) (1l| , R2 = |0l) (- + +1 + (+ - -I , 



i?3 = 10, 



+ I + IU) 



, R4 = \0i 



II, 



(30) 



Using simple algebra, it turns out that the 8x8 matrix representation [Ri] with I = 1,..,4 of the recovery operators 
is given by. 



[Ri] = Ell + Ess, [R2] — E12 + Est, [R^] = i?i3 + Esq, [i?4] = -f- Esz, 



(31) 



where Eij is the 8x8 matrix where the only non-vanishing element is the one located in the ij-position and it equals 
1. The action of this recovery operation TZ on the map A'^'^^ (p) in (fTO|) leads to. 
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ALl™, ip) ^ 7^oA(3) (p) ^ ^ (i?, a;,) p (i?, A',) 



(32) 



fe=0 /=1 



We want to describe the action of TJ.oA^'^^ restricted to the code subspace C. Therefore, we compute the 2x2 matrix 
representation [i?/j4y|^ of each RiA'^ with I = 1,.., 4 and k = 0,.., 7 where, 



[RiMc ~ { {1l\RiA',\0l) {1l\RiA',\W 



(33) 



Substituting (|26p and pop into p3p . it turns out that the only matrices with non- vanishing trace are given 

by. 



[RiA'o_ 



J'o I 1 j ' ^^2^ij|c - yPi I 1 



^3 I 1 



(34) 



Therefore, the entanglement fidehty ^p^iso p) defined as, 
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-^phL (a^. p) - (2/ - + p) + ^ (-4/ + - 2p) + (2/ _ Sp^ + l) . (36) 



results. 



DFS for correlated phase flips. By considering the case of (|16p with n = 2 qubits it follows that the error superop- 
erator A associated to channel is defined in terms of the following error operators, 

A ^ {A'„.., A',} with A(2)(p) E ^^-K^ ^^d, E ^'J^'fc = ^4x4. (37) 

k=0 k=0 

In an explicit way, the error operators {A'q,-., A'^} are given by, 

A'o = VpF^' ® ^'i = ^/W^z^ ® /2, = y^/i ® = ^/^z^ ® z^, (38) 

where the coefficients p^^-* for fc = 0,.., 3 read, 

-(2) ~(2) ~(2) -(2) .„„s 

Po = POOPO, Pi = PlOPO, P2 = PQlPl, Ps = PllPl- (39) 

We encode our logical qubit with a simple decoherenee free subspace of two qubits given by , 

|0) \0l) |01) and, |1) |U) |10) . (40) 

The set of error operators satisfying the detectability condition ^o|, PcA'f^Pc ~ ^A' Pc, where Pq — (Ol| + 
is the projector operator on the code subspace C = Span{\OL) , is given by, 

•^detectable 

{A',, A',} c A. (41) 

Furthermore, since all the detectable errors are invertible, the set of correctable errors is such that 
-^Lrrectabie-^correctabie IS detectable. It foUows then that, 

-^correctable — -^detectable 

C A. (42) 
The action of the correctable error operators -4correctabie on the codewords |0l) and jl^) is given by, 

|0i) A, |0i) = VpFiOI) , A', |0i) = -VpFiOI) , \H) ^ A', \W = VpFiiO) , A', \U) = -y^f |10) . (43) 

The two one-dimensional orthogonal subspaces V"^ and V^^ of "H^ generated by the action of -4corrcctabie on |0l) 
and are given by, 

V"^ = Span{\vl'') = |01)} and, V^^ = Span{\vl'-) = |10)} . (44) 

Notice that ® V^^ ^ This means that the trace preserving recovery superoperator TZ is defined in terms of 

1 

one standard recovery operator Ri and by the projector R±^ onto the orthogonal complement of ^ V*^ , i. e. the 

part of the Hilbert space which is not reached by acting on the code C with the correctable error operators. In 
the case under consideration, 

2 

def I 



i?i |01) (Oil + |10) (10| , i?± = E Vs) {rs\ , (45) 

where is an orthonormal basis for (V°^ ® V^'^)''". A suitable basis '^^yo^.^yii,')^ is given by, 

^(vOr.^v^r.y={^i = m,r, = \U)}. (46) 
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FIG. 1: Threshold curves for code efTectiveness: concatenated code (dashed hne), DFS (sohd hne). The concatenated code 
works in the parametric region below the dashed line. The DFS works in the region above the solid line. The repetition code 
works for any value of ^ when p < 0.5. 



The action of this recovery operation TZ vi^itli i?2 = R± on the map A*^^^ {p) in ([57]) yields, 

Ag)_, (p) ^ (7^oA(2)) (p) ^ ^ {R,Ai) p {RiA',)^ . (47) 



fc=0 1 = 1 



We want to describe the action of T^oA'^^ restricted to the code subspace C. Therefore, we compute the 2x2 matrix 
representation of each RiA'^. with ^ 1, 2 and k = 0,.., 3 where, 

. clef / (OL|i?,AyOL) {Ql\RiA',\Il)\ 

Substituting ([38)) and pS]) into p8)) . it turns out that the only matrices |^ with non- vanishing trace are given 



by. 



[^i^olic = ( J ? ) ' [^i^3]|c = - V5F ( J ? ) (49) 



(2.) _____ 

Therefore, the entanglement fidelity J-jjpg {p, p) defined in (|10p results, 

t'^Is {p, p)^p (-2p2 + 2p) + (2p2 - 2p + 1) . (50) 

We point out that error correction schemes improve the transmission accuracy only if the failure probability V {p, p) 
is strictly less than the error probability p [32| , 

rip,p) = i-T{p,p)<p. (51) 

The threshold curves for code effectiveness p (p) arc defined by the relation 1 — T {p (p) , p) ^ 0. They allow to select 
the two-dimensional parametric region where error correction schemes are useful. The threshold curve for the DFS 
for the model appears in Figure [T] The DFS works only in the parametric region above the thin solid line in Figure 
[U while the repetition code works for all values of the memory degree p when the error probability p is less than 0.5. 

Concatenated Code. Following the line of reasoning presented in |33| . we consider the concatenation of the two 
codes just described. In the case of (fTB]) with n = 6 qubits and correlated errors in a dephasing quantum channel, 

1 

A^^\p) — Pig\i^Pi^\i^Pi^\i^Pi,j\i2Pi2\iiPhiAigAi^Ai^Ai^Ai.-^Ai-^)p{AigAi^Ai^Ai^A , (52) 

ii,..., 16=0 
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The error superoperator A associated to channel (|52p is defined in terms of the following error operators, 

A ^ {A'o,.., A',,} with A(6)(p) A'kP^l and, ^ = /64X64. (53) 

fc=0 k=0 

The error operators in the Kraus decomposition ((53)) are X]fc=o (D = 2^ = 64, where (^) is the cardinality of weight-fc 
error operators. 

We encode our logical qubit with a concatenated subspace obtained by combining the decoherence free subspace in 
(|40)) (inner code, Cdfs = Cinncr ) with the repetition code in ([23)) (outer code, Cphase = Couter)- We obtain that the 
codewords of the concatenated code C = Cdfs ° Cphaso arc given by, 

|0l> = 1+ + + - — > , lU) I + ++> . (54) 

Recall that the dctectability condition is given by PcA'f.Pc — Xa' Pc where the projector operator on the code space 
C is Pc = |0l) (Ol| + |1l) (1l|. Observe that, 

PcA'^Pc = {i)MQL) |0l) (OlI + {QMW |0l) (IlI + IU) (Oil + {Il\A',\\l) |U) (UI • (55) 

Therefore, it turns out that for detectable error operators we must have, 

{QL\A'k\QL) = (IlI^I-IIl) and, {Ql\A'^\Il) = {Il\A'^\Ql) = 0. (56) 

In the case under consideration, it follows that the only error operator (omitting for the sake of simplicity the proper 
error amplitudes) not fulfilling the above conditions is proportional to, 

Z^Z^Z^Z^Z^Z"^. [hi) 

For such operator, we get 

{^l\Z^Z''-Z'^Z^Z^Z^''\\l) = 1 7^ and, (\i\Z^ Z"^ Zi"^ Z^ Z'^ Z^^^d^) = 1 7^ . (58) 

Therefore Z^Z^Z^Z^Z^Z^ is not detectable. Thus, the cardinality of the set of detectable errors v4dotoctabie is 63. 
Furthermore, recall that the set of correctable errors -^correctable SUCll tllB-t -^.^^^.^^^^.^j^i^vAcorrectable detectable (in 
the hypothesis of invertible error operators). Therefore, after some reasoning, we conclude that the set of correctable 
errors is composed by 32 error operators. The correctable weight-0, 1 and 2 correctable error operators are (omitting 
the proper error amplitudes), 

{nweight-o . {Z\ Z\ Z\ Z^)^^^^^^_^ , (59) 



and, 

jz^z^, z^z^, z^z^, z^z^, z^z^. z'^z^, z'^z^, z:^z^. z''-z^. Z^Z^, Z^Z^, Z^Z^, Z'^Z^. Z^Z^, Z^Z^] . ^ , 

L '''''''''''''' J weignt-2 ' 

(60) 

respectively. The correctable weight-3 errors are, 

{Z^Z^Z^, Z^Z^Z^\ Z^Z^Z^, Z^Z^Z^, Z^Z^Z^, Z'^Z'^Z'', Z^Z^Z^, Z'^Z^Z^, Z^Z'^Z^, Z'^Z^Z^}^^.^^^_^ . (61) 

There arc no wcight-4, 5 and 6 correctable error operators. The action of the correctable errors on the codewords in 
([S3]) is such that the Hilbert space H2 can be decomposed in two 32-dimensional orthogonal subspaces and V^^. 
In other words, "Hj = © ^'^^ where 

V"- = Span \ \vl-^,) = ^L=A', \0l) \ and, V'- - Span \ \vll,) = -^A', |U) I , (62) 
I J I J 

with Af^ e .4coiTcctabic Vfc = 0,..., 31 (numbering the correctable error operators from to 31). Notice that (w^^ \vj,f \ — 



ij. ^ ^correctable — "v, i^iiuiiiuvji iiig uiic i.vjiicv^tciuic ciiui vjpcia,t<jj.o liuiii u uuux;. i^ijun.,c tiicit yi/j, | i/j,, 

Skk'Sij, with k, k' £ {0,..., 31} and i, j £ {0, 1} since, 

\ yPk yPk' I yPk pV vpI Pi' 



where we have used the fact that the square (Hermitian) matrix a'^^, equals y p\, p\., Skk' ■ The recovery superoperator 
TZ o {Ri} with I = 1,., 32 is defined as [IH, 
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-^Ehr)(«ri, (64) 

1=0 

where the unitary operator Vi is such that Vi = for i G {0, 1}. Notice that, 

1 

Ri Vi J2 hi") ("^ I = |0l> ("?" I + |1l> I • (65) 



It turns out that the 32 recovery operators are given by, 



Ri+i = i?i4if = (|0l) (OlI + (IlI) (66) 



'Pi VPi 
with / G {0,..., 31}. Finally, the action of this recovery operation TZ on the map A*-^^ (p) in ([53]) leads to, 

2*^-1 32 

MlLer (P) ^ (7^oA(6)) (p) E (^'^'^O P (^'^^^)^ • (67) 



fe=0 1 = 1 



We want to describe the action of TJ.oA'^-' restricted to the code subspace C. Recalling that A'^ = A[\ it turns out 
that, 

{zl\Ri+iA',\jl) - ^ (*L|Oi) (Ol\4a[\jl) + (*l|1l) (uIAU'.Ijl) . (68) 

We now need to compute the 2x2 matrix representation of each RiA!^ with / = 0,.., 31 and k ~ 0,.., 2^ — 1 
where, 

. dcf / (OL|i?,+iA',|OL) {^l\Ri+iA!^\\l)\ ,.qx 

[^/+i^.-J|c - (u|i?,+,A'jOi) (U|i?,+iAyu) ; ■ ^^^^ 

For I, k = 0,.., 31, we note that becomes, 

[o,\4a',\o,) \ r;. fi 

while for any pair {I, k) with I — 0,.., 31 and k > 31, it follows that, 

{OlIRi+iA'M + {lL\Ri+iAi\lL) = 0. (71) 
We conclude that the only matrices with non- vanishing trace are given by [i?;+i^J]|^ with I = 0,.., 31 where, 

J). (72) 

Therefore, the entanglement fidelity J-^onc (m> p) defined in becomes, 

•^ionc (a*: p) = PooPo + 2pooPioPo + ^ploPoiPioPo + pIqPioPiiPo + 3pIoPoipIoPo+ 

+ P00P01P1QP1 + iploPoiPioPnPo + ipooPoiPwPo + ^pIqpIiPioPi + pIoPoiPiiPi + 

+ PoiPiaPQ + 2pooPoiPioPi + PaoPaiploPiiPo + PoaPoiPwPiiPi + PqqPqiPiqpIiPq+ 

+ '^PoiPwPiiPo + PooPoiPioPiiPi + PoiPioPi- (73) 
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FIG. 2: Threshold curves for code performance: concatenated code vs. DFS (thin sohd hne) and concatenated code vs. 
repetition code (dashed Hne). The concatenated code outperforms the DFS in the parametric region below the solid line and 
outperforms the repetition code in the region below the dashed line. 



Substituting p2|) into (O, we finally get 

•^ioL (/^, P) = (- V + 11/ - lOp' + 3p2) + (10/ - 25/ + 22/ - 8p^ + p) + 

+fi^ (-6/ + 12p3 - 7p2 + p) + ^2 (^_20p'5 + 58/ - 60p3 + 25^^ - 3p) + 

+fi {20p^ - 53/ + 46p3 - 13p2) + (-6p^ + 15/ - lOp^ + l) . (74) 

The threshold curve for code efFectivencss concerning the concatenated code defined in (j54p for our noise model 
appears in Figure [TJ It turns out that the concatenated code vi^orks in the parametric region below the dashed line. 

To uncover the parametric region where one code (say, code-1) outperforms another code (say, code-2), we consider 
the threshold curves for code performances /2 (p) defined by the relation ^code-i {p- (p) , p) — -7^;odc-2 {p- (p) , p) = 0. 
We emphasize that in view of equations ((36)) . (|50p and ((74)) . it turns out that the concatenated code outperforms 
the DFS in the parametric region below the solid line (J^ionc {p- ip) , p) — J^dfs (P- (p) ; p) = 0) and outperforms the 
three-qubit repetition code in the region below the dashed line (J^conc {p (p) , p) — -T^phLc iP (p) j P) = 0) Figure [2] 
For the sake of clarity, in Figure [3] we plot the entanglement fidelities ((36)) (thin solid line), ([74]) (dashed line) and 
([SH]) (thick solid line) for p = 10~^. Our analysis explicitly shows that none of the two codes (DFS and repetition 
code) is effective in the extreme limit when the other is, the repetition code still works for correlated errors, whereas 
the error avoiding code does not work in the absence of correlations. Finally, our final finding leads to conclude that 
there is a parametric region characterized by intermediate values of the memory parameter where the concatenated 
code in ([M)) is particularly advantageous (see Figure 



B. Code entropy-based analysis 

Repetition code for correlated phase flips. The error correction matrix Fphase for the three-qubit repetition code 
considered is given by, 

/pIoPo 

J. _ pooPwPo 
P'^"'^'^ ~ PoiPwPo 

V PooPoiPi 
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FIG. 3: Entanglement fidelity vs. memory parameter fj, for p = 10 ^: concatenated code (dashed line), DFS (thick solid line) 
and repetition code (thin solid line). 



Therefore the repetition code entropy 5rc (/^, p) results from (fT4|). 

^Rc = -pIoPo log2 (PooPo) - PooPioPo log2 {pooPwPq) - PaiPiQPa logs (poiPioPo) - PooPoiPi log2 {paoPaiPi) , (76) 



with. 



PooPo 



P 



(77) 



In this case it turns out that there is no value of the memory parameter p for which Crc is a unitarily correctable 
code. In other words. Sue (m, p) for any < /i < 1 with < p < 1. Therefore, it never occurs the case 
where the effort for recovering the code is minimum. For instance, in the extreme limit of /i = 1 it follows that 
iSrc (a* = p) ~ ^ (1 ^ p) log2 (1 ^ P) 7^ 0- From Figure 21 we notice that in general in the limit of small error 
probabilities (p <C 1), the nearness of Cflc to a DFS increases when p increases and/or the error probability p 
decreases. Therefore the RC entropy analysis, although not particularly enlightening, confirms that there is no pair 
of parametric values /.t and p for which Crc provides a fully protected space in the presence of correlated phase-flip 
errors and the effort required for recovering the quantum state corrupted by the phase-flip noise increases when p 
increases. 

DFS for correlated phase flips. The error correction matrix Fdfs for the DFS considered is given by. 



where. 



DFS 



dcf , dcf dcf ; 

a = PooPo, = piipi, c = -^JpmPoPiiPi- 



Diagonalizing Adfsi it turns out that. 



with. 



[AdFs] diagonal ^ ( Q X_ 



A+ = ^ 



[a + b)± ^J[a- bf + 4c2 



(78) 



(79) 



(80) 



(81) 
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FIG. 4: RC entropy vs. /i for p = 2 x 10 ^ (thick solid line), p = 2.5 x 10 ^ (thin solid line) and p = 3 x 10 ^ (dash line). 
The RC entropy does not vanish for /j. = 1. 

Substituting ([79l) into ([8T|) and recalling that, 

Po = (1 - P) , Pi = P, Poo = ((1 - m) (1 - p) + m) , 

Poi = (1 - m) (1 -P) , Pio = (1 - m)p, Pii = ((1 - /i)p + /^) , (82) 

we obtain, 

A+ = (-2p^ + 2p) + (2p2 - 2p + 1) and, A_ = 0. (83) 
In conclusion the code entropy 5dfs (p-, p) results from (|14p. 

5dfs (m, P) = - [m (-2p' + 2p) + (2p2 _ 2p + 1)] log2 [fi {-2p^ + 2p) + {2p^ - 2p + l)] . (84) 

From Figure HI it follows that Cdfs is a unitarily correctable code only in the extreme limit of /i = 1. In such a 
case, SuFS (m = 1; p) = for any < p < 1 and thus the effort for recovering the code is minimum since a trivial 
identity recovery operation will do the job. From Figures [5] and [6l we notice that in general in the limit of small error 
probabilities (p <SC 1), the effort required for recovering the quantum state corrupted by the phase-flip noise increases 
when /i increases and/or p increases. Thus, the DFS entropy analysis leads to the conclusion that Cdfs provides an 
especially useful error correction scheme for very small error probabilities and highly correlated phase-flip errors. 

Concatenated Code. After some straightforward but tedious algebra, it turns out that the entropy of the concate- 
nated code iSconc. i^J', p) is given by. 



31 

5conc. (m, p) = - X! -^J ^°S2 f] (m, P) , 

J=0 



(85) 
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FIG. 5: DPS entropy vs. p for j-i = 1 (thick solid line), /i = 0.90 (medium solid line), jj, = 0.75 (thin solid line), fi = 0.5 (thick 
dash line) and fi — (thin dash line). 




FIG. 6: DFS entropy vs. memory parameter for p = 10 ^ (dash solid line), p = 5 x 10 ^ (thin solid line) and p = 10 
(medium solid line). The DFS entropy vanishes for fi = 1. 
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FIG. 7: Concatenated code entropy vs. memory parameter for p — 10~^ (dashed solid line), p = 5 x 10~^ (thin solid line), 
p = 10~^ (medium solid line), 

where, 

fo = PqqPo, /i = /2 = PooPioPo, fs = fi = h = fe = pIqPoiPioPo, fr PooPioPiiPo, 
/s = /g = /lo = PooPoiP?oPo, /ii = PooPoiPioPi, /12 = /13 = /u = PooPoiPioPiiPo, 

2 2 2 2 3 

.fl5 = /16 = ,/l7 = POOPoiPloPO, ./18 = /l9 = ./20 = PooPoiPlOPl, /2I = PooPolPllPl, 

/22 = PoiPioPo, /23 = /24 = PoopIipIqPi, hb = PooPoiPwPilPo, he = ploPoiPioPiiPi, 

/27 = PooPoiPioPnPo, /28 = /29 = PqiPioPhPo, /30 = PoopliPioPiiPi, hi =pIiPwPi. (86) 

and, 

Pa = (1 - P) , Pi = Poo ((1 - /") (1 - P) + m) , 

Poi = (1 - m) (1 -p) , Pio = (1 - Pu = ((1 - P)P + P) ■ (87) 

It can be easily checked that in the extreme limit of = 1 it follows that 5conc. (m = 1j P) = ^ (1 ^ P) log2 (1 ^ P) 7^ 0. 
In general, it appears that results similar to those obtained within the RC code entropy analysis hold (see Figure [7]). 
Then, it seems that the nearness of Cconc. to a DFS increases when ^ increases and/or the error probability p decreases. 
In particular, the effort required for recovering the quantum state corrupted by the phase-flip noise increases when p 
increases. 

Our analysis for Model I allows us to conclude that the entanglement fidelity is generally a better code performance 
quantifier than the code entropy since it is easier to handle computationally and, most of all, it allows to compare 
performances of different error correction techniques applied to the same error model in a smoother way. However, 
it appears that the code entropy is especially informative when quantifying the performance of DFSs in the limit of 
highly correlated and very low error probabilities. 

IV. MODEL II: ASYMMETRIC DEPOLARIZING NOISY QUANTUM MEMORY CHANNEL 

In this Section, wc consider both symmetric and asymmetric depolarizing noisy quantum memory channels and 
perform quantum error correction via the five qubit stabilizer code C[[5 1^3]]. We characterize this code by means of 
the code entropy and the entanglement fidelity as function of the error probability and the degree of memory. In 
particular, we uncover that while asymmetric depolarizing errors do not affect the entanglement fidelity of the five 
qubit code, they do affect its code entropy. 
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A. Entanglement fidelity-based analysis 

For the symmetric case, consider five qubits and Markov correlated errors in a depolarizing quantum channel A*-^-* (p), 

3 

Pi^liiPiil isPia \i2Pi2\i1Pi1 ^'h ^'h P^l^ ^Ig ^Is 

il, 42, is, ^4, 25—0 

where Aq = I, Ai = X, A2 =Y, A3 = Z are the Pauli operators. The coefficients (conditional probabilities) 

with I, m ^ {0, 1,..., 5} satisfy the normalization condition, 

3 

E P^r,\iiP^i\^3P^3\'i.2P^2\^lPil = 1- (89) 

ii, 12, is, U, 15=0 

For the depolarizing channel k^^\p), coefficients are considered as, 

Pk\k-i'=^ - P')Pk + ^^5k\k~l, Pfe=o = 1 -P, Pfc=i, 2, 3 =p/3. (90) 
Following [ll[, it turns out that the entanglement fidelity -T^iymmeLic (/^' P) results from pUj) 

15 

•^iymmLc (M, P) = E (^' = /o + (A + - + U) + (/t + - + As) , (91) 

fe=0 

where the functions (^, p) are given by, 

/o = PooPo, /i = ••• = /e = PooPioPo, /7 = ••• = /i5 = PooPoiPioPo, (92) 

with, 

P 

PO = 1 - P, Pi = P2 = P3 = 2 , Poo = (1 - ^) (1 - P) + A^, 

POI = P02 = P03 = (1 - (1 - P) , PlO = P20 = P30 = 2 (1 - A^) • (93) 
For the asymmetric case, we assume that the error probability p may be written as, 

P = Px +PY +Pz, (94) 

where, 

Px = axp, Py = ayp, pz = azp, (95) 
with ax + tty + ctz = 1. Notice that in the symmetric case, we simply have ax = cty = az = -I- It turns out that 
the -FiLyifiitric (M- p) is given by [HI, 

15 



•^l„ctric iP^ P; ax, ay, az) ^Y.^k (/^' = /o + (/( + - + /e) + (/y + - + fi^) > (96) 



-[[5,1,3]] 

;ric \h^^ ^"1 "A , «y , "z; — 

fc=0 

where the functions (^, p; ax, ay, az) read, 

/o = PooPo; f'l = PooPoPio, /2 = P00P0P20, /s = P00P0P30, f'i = P00P01P17 /s P00P01P2, /e P00P01P3, 

/? = /s = /o = PooPoiPoPio, fw = fii = .fi2 = P00P01P0P20, /(s = fiA = lib = P00P01P0P30, (97) 

with, 

Po = 1 - P^ Pi = "XP, P2 = Olyp, P3 = azP, Poo = (1 - M) (1 - P) + 
POI = P02 = P03 = (1 - Ai) (1 - P) , PlO = axp (1 - m) > P20 = Oiyp {I - p) , P30 = "zP (1 - /") ■ (98) 
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FIG. 8: Threshold curve for the five-qubit stabihzer code effectiveness. The C[[5^i^3]] code works in the parametric region below 
the curve. 
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FIG. 9: 

vs. n with < fi < 0.33 ( for fi > 0.33, the error correction scheme is not effective anymore) for 
p = 4.50 X 10"^ (thick solid line), p = 4.10 x 10~^ (thin solid line) and p = 3.75 x 10"^ (dashed line). 

Recalling that in the symmetric case Pi — P2 — Pa ~ ^ and pio = P20 ~ P30 — § (1 ^ m) and substituting (|98p in (|96p . 
it follows that 

•^ilymmctric (/^' P'^ Otx^ay, az) = ^^y^^ln^ (M, P) ■ (99) 

Therefore, we conclude that the performance of the five-qubit stabilizer code, when quantified by the entanglement 
fidelity, remains unaffected by the asymmetry of the depolarizing error probabilities. From (|96|) it can be shown that 
such error correction scheme only works for low values of /i (see Figure [8]). Furthermore, it also turns out that the 
performance of the five qubit quantum stabilizer code is lowered by increasing values of the degree of memory ^ and 
increasing values of the error probability values p (see Figure [5]). For a more detailed overview of such findings, we 
refer to 
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FIG. 10: DifTerence between the five qubit stabilizer code entropies <Sjyj-|jnf|]|ri;. (p) and 5^[^'^'^il 
tions (/i = 0) with ~ ciy = j and = i. 



asymmetric 



(p) in the absence of correla- 



B. Code entropy-based analysis 

Omitting technical details and following the work presented in the previous Section, it turns out that the five qubit 
code entropies for the symmetric and asymmetric depolarizing quantum memory channels are given by, 



15 

5!!!:ifiL,„ (^, p) = - ^ fj (m, p) loga fj (m, p) 



symmetric ' 



-0 



(100) 



and, 



5[[5,1,3]] 
asymmetric 
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(/i, p; ax, ay, az) = - ^fj (M: P; <^x, ay, az) logs /j (M: P\ "x, ay, az) 



(101) 



respectively. The explicit expressions for the functions (/i, p) and /j (/i. p; ax, ay, az) are given in and (|97p . 
respectively. From (llOOp and (jlOip . it can be shown that 



'^iymm^ilic (M, P) ^ 4Lymmotric (M, P; , ay , az) ■ 



(102) 



Numerical evidence of the disequality in (|102p can be easily verified at least for a suitable choice of model parameters 
/.t, ax, ay and az- In Figure [TUl we plot the change of the five qubit code entropy A5 (p) as function of the 
error probability p. Specifically, A5 (p) is the difference between the five qubit code entropies ^jlyj^j^j^]!^;^. {p) and 



S. 



[[5,1,3]] 
asymmetric 



(p) in the absence of correlations (/x = 0) with ax = ay 



\ and az = \- The positivity of A5 (p) leads 



to conclude that when the error model considered is defined by a suitable choice of numerical values of the model 
parameters, it can happen that the effort required for recovering the five qubit code in the asymmetric case is less than 
in the symmetric case. This scenario never occurs when quantifying the performance of the five qubit code applied 
to asymmetric and correlated depolarizing errors by means of the entanglement fidelity. 

From an analytical point of view, (|102p can be readily understood by noticing that {[i, p) in (|100p and 

{fi, p; ax, ay, az) in (|101|) are sums of terms that are nonlinear in the functions (amplitude square 



[[5,1,3]] 
asymmetric 



5, 

probabilities of the correctable error operators) fj (/i, p) and /j (/x, p; ax, ay, az), respectively. As a matter of fact 
it turns out that. 



15 



15 



(103) 



k=0 



k=0 
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with f^ = fo= ptoPo, f[ + ... + = /i + ... + /6, + ... + /Is = /7 + ... + /i5 and, // ^ /j for i, j = 1,..., 6; 
fr = fa = fg ^ PoaPoiPoPiol fw = f'li = /12 = and /(g = = /(g = PooPoiPoPao. However, the 

constraint in (|103p together with the nonUnearity in / and /' characterizing the expressions of the code entropies in 
(ITOO)) and ((TOTI) imply that, 

15 15 

E/fcl°g2/fc^E/fcl0g2/^ (104) 

fc=0 fe=0 

In conclusion, since the entanglement fidelity is a linear combination of the amplitude square probabilities of the 
correctable error operators, Eq. (|103p holds true, and the equality in ([QQ]) is proven. On the contrary, the code 
entropy is a nonlinear combination of the amplitude square probabilities and the use of (|104p leads naturally to the 
inequality in (|102p . 

Therefore, we conclude that the performance of the five-qubit stabilizer code, when quantified by the code entropy 
and not by the entanglement fidelity, does remain affected by the asymmetry of the depolarizing error probabilities. 
Our finding provides a neat quantitative manifestation of the conceptual fact [26j that no single code performance 
quantifier holds all the information on a code. 



V. FINAL REMARKS 



In this article, we studied the properties of error correcting codes for noise models in the presence of asymmetries 
and/or correlations by means of the entanglement fidelity and the code entropy. We considered a dephasing Markovian 
memory channel (Model I) and both symmetric and asymmetric depolarizing quantum memory channels (Model II). 
For each model, we presented both an entanglement fidelity-based and a code entropy-based analyses. For Model I, 
we used three codes: the repetition code C_rc, the DFS Cdfs and the concatenated code Cdfs ° Crc- For Model II, 
we employed the five qubit stabilizer code C[[5^i.3]]. 

For Model I, the entanglement fidelity-based analysis allows to find out the parametric regions where the chosen 
error correction schemes are effective (see Figure [T]). This analysis is also suitable to determine where, within such 
parametric regions, one code outperforms the other (see Figure [5]). In particular, we showed that the concatenated 
code quantified by the entanglement fidelity is useful to combat partially correlated phase-flip errors (sec Figure 
|3]). The code entropy-based analysis for Model I leads to the conclusion that the effort required for recovering the 
quantum state corrupted by correlated phase-fiip errors increases when the memory parameter p, decreases and the 
error probability p increases (see Figures 2] and [5]) . Furthermore, it turns out that only Cdfs is a unitarily correctable 
code for the noise model considered in the limiting case of /i = 1 (see Figure [5]). When applied to the concatenated 
code, the code entropy-based analysis is not as enlightening as the entanglement-fidelity based analysis. It is not 
particularly useful for drawing performance comparisons with the repetition code (see Figure [7]). 

For Model II, the entanglement fidelity-based analysis allows us to find out that the five qubit stabilizer code only 
works for small values of fi (see Figure |8]). Furthermore, the performance of this code applied to symmetric and 
asymmetric correlated depolarizing errors is lowered by increasing values of the memory parameter fj, and increasing 
values of the error probability p (see Figure ^ . Finally, the code entropy-based analysis applied to Model II leads to 
an interesting result. While asymmetry in the depolarizing errors does not affect in any case the performance of the 
five qubit stabilizer code quantified by means of the entanglement fidelity, it may affect positively the performance of 
C[[5,i,3]] by lowering the effort required for recovering the code subjected to asymmetric depolarizing errors (see Figure 

Although our work is limited to only few error models in the presence of correlations and asymmetry and we 
only perform error correction by means of few quantum codes, we feel we have gathered enough evidence to draw 
the following conclusions on the code performance quantifiers employed: a) while the code entropy may capture 
new undetected features of a code with respect to those encoded into the entanglement fidelity and deserves further 
investigations, it is certainly harder to compute for realistic physical models and non-trivial quantum codes; b) the 
code entropy seems to lack the reliability and practicality characterizing the entanglement fidelity. 

In summary, the entropy of a code appears to be an interesting auxiliary tool for quantifying the code performances. 
In particular, it is our intention to explore in future investigations the possibility of devising a hybrid code performance 
quantifier whose nonlinear structure in the amplitude square probabilities of the correctable error operators may not be 
necessarily characterized by the logarithmic behavior which identifies the code entropy. For the time being, motivated 
by our analysis and in agreement with j2^ [s^ . [s^ , we believe that the entanglement fidelity remains the most relevant 
tool to maximize in schemes for quantum error correction. 
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